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Strong Disorder Renormalization for the Random Transverse Field Ising model leads to a com- 
plicated topology of surviving clusters as soon as d > 1. Even if one starts from a Cayley tree, the 
network of surviving renormalized clusters will contain loops, so that no analytical solution can be 
obtained. Here we introduce a modified procedure called 'Boundary Strong Disorder Renormaliza- 
04 ■ tion' that preserves the tree structure, so that one can write simple recursions with respect to the 

number of generations. We first show that this modified procedure allows to recover exactly most 
of the critical exponents for the one-dimensional chain. After this important check, we study the 
RG equations for the quantum Ising model on a Cayley tree with a uniform ferromagnetic coupling 
J and random transverse fields with support [hmin, hmax]- We find the following picture (i) for 
^ J > hmax, only bonds are decimated, so that the whole tree is a quantum ferromagnetic cluster (ii) 

. for J < hmin, Only sites are decimated, so that no quantum ferromagnetic cluster is formed, and the 

QQ ' ferromagnetic coupling to the boundary coincides with the partition function of a Directed Polymer 

model in a random medium (iii) for hmin < J < hmax , both sites and bonds can be decimated : the 
quantum ferromagnetic clusters can either remain finite (the physics is then similar to (ii), with a 
quantitative mapping to a modified Directed Polymer model) or an infinite quantum ferromagnetic 
Ch ' cluster appears. We find that the quantum transition can be of two types : (a) either the quantum 

C , transition takes place in the region where quantum ferromagnetic clusters remain finite, and the 

ly-j ■ singularity of the ferromagnetic coupling to the boundary involves the typical correlation length ex- 

ponent i/typ = 1 (b) or the quantum transition takes place at the point where an extensive quantum 
ferromagnetic cluster appears, with a correlation length exponent v ~ 0.75. 
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INTRODUCTION 



^ ■ In this paper, we consider the quantum Ising model defined in terms of Pauli matrices 

o 

^- n^-Y.J.,,a:a^-Y,h,a: (1) 

<i,j> i 

(N 

^ where the transverse fields hi > are independent random variables drawn with some distribution 7r(/ii) and where 
the nearest-neighbor couplings take a non-random value Ji,j = J > (of course, one can also consider the case where 
the couplings are also random, but we have chosen to restrict to the case of non-random J to simplify a bit the 
notations). In dimension d — 1, exact results for a large number of observables have been obtained by Daniel Fisher 
[l[ via the asymptotically exact strong disorder renormalization procedure (for a review, see 0) : the transition is 
1^ ' governed by an Infinite-Disorder fixed point and presents unconventional scaling laws with respect to the pure case. 

In dimension d > 1, the strong disorder renormalization procedure can still be defined. However, it cannot be solved 
analytically, because the topology of the lattice changes upon renormalization, but it has been studied numerically 
with the conclusion that the transition is also gove rned by an Infinite-Disorder fixed point in dimensions d = 2,3,4 
[sl-fisj and on Erdos-Renyi random graphs l3- These numerical renormalization results are in agreement with 



^> ' the results of independent quantum Monte-Carlo in d = 2 [IJ, [15| . 
5^ , Nevertheless, the complicated topology that emerges between renormalized degrees of freedom in dimension d > 1 
tends to obscure the physics, because a large number of very weak bonds a re g enerated during the RG, that will 



eventually not be important for the forthcoming RG steps. In a recent work 16|, we have thus proposed to include 
strong disorder RG ideas within the more traditional fixed-length-scale real space RG framework that preserves the 
topology upon renormalization : the idea is to allow some changes in the order of decimations with respect to the full 
procedure in order to maintain a simple spatial renormalized lattice. We have checked that in dimension d = 1, this 
procedure correctly captures all critical exponents except for the magnetic exponent /3 which is related to persistence 
properties of the full RG flow. In the present paper, we introduce a similar idea for the case of the Cayley tree 
geometry. If one applies the full strong disorder RG to a Cayley tree, the tree structure is rapidly destroyed and the 
network of surviving clusters contains loops. Here we thus introduce a modified procedure called Boundary Strong 
Disorder Renormalization that preserves the tree structure, so that one can write simple recursions with respect to 
the number of generations. The 'price' is again that the magnetization will not be well taken into account, but one 
can hope that other scalings are correctly captured, as in fl^. We will show that this is indeed the case in dimension 
d = 1 by a direct comparison with the exact solution of the full RG procedure. For the Cayley tree, we find that when 
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only sites are decimated within the Boundary Strong Disorder RG, one obtains a quantitative mapping towards some 
Directed Polymer on the Cayley tree : this relation with the Directed Polymer has been already obtained for the Cayley 
tree via some approximations within the Quantum Cavity Approach [iTHlSj, and for arbitrary networks via simple 
perturbation deep in the disordered phase l2C| . However here we also consider the possibility of bond-decimations to 
build quantum ferromagnetic clusters which may become important near the transition. 

The paper is organized as follows. In section [III we introduce the Boundary Strong Disorder Renormalization 
procedure for a Cayley tree of branching ratio K (i.e. each non-boundary site has {K -t- 1) neighbors). In section Hill 
we show that for the one-dimensional chain corresponding io K = 1, this procedure is able to reproduce the most 
important critical exponents. After this important check, we study the RG equation for real trees having K > 1, as 
a function of the value of the ferromagnetic J with respect to the support [hmimhmax] of the distribution 7r(/ii) of 
random fields : the case J > h^ax where only bonds are decimated, the case J < h^in where only sites are decimated, 
and the case hmin < J < h„iax where both bonds and sites can be decimated are discussed in sections IIVI fVl and IVII 
respectively. 

II. BOUNDARY STRONG DISORDER RENORMALIZATION PROCEDURE 

As recalled in Appendix [Xj the Strong Disorder Renormalization for the quantum Ising model of Eq. [1] is an 
energy-based RG, where the strongest ferromagnetic bond or the strongest transverse field is iteratively eliminated. 
In this section, we introduced a modified procedure, called Boundary Strong Disorder RG, that preserves the tree 
structure, in order to write explicit recursions for probability distributions of renormalized variables. 

A. Notations 




ext 



FIG. 1: Renormalized structure at RG step n for a Cayley tree of branching ratio K — 2. All sites of generations m > n + 1 
have not yet been modified with respect to the initial model, i.e. the sites are characterized by their initial random fields 
/im'(z), and are connected to their neighbors by the initial ferromagnetic coupling J. All sites of generations with m < n — 1 
have disappeared. All sites i of generation n have renormalized transverse fields h^{i) and are connected to the formal external 
spin (called "ext") via some renormalized coupling J^'^^^{i). 

We consider a Cayley tree of branching ratio K with N generations : we label with n = the generation of external 
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leaves that are formally connected to an exterior spin with the coupling J^^^q-^ , in order to keep track of the coupling 
to the boundary upon decimation. The next generations are labeled with n = 1, 2,3, .. that measures the distance to 
the boundary. 

The RG procedure consists in the renormalization of generations in the order of the label n starting from n = 0. 
At RG step n, we have thus the following structure (see Fig. [T]) 

- all sites of generations ni > n + 1 have not yet been modified with respect to the initial model, i.e. the sites are 
characterized by their initial random fields h™^{i), and are connected to their neighbor by the initial ferromagnetic 
coupling J. 

- all sites of generations with m < n — 1 have disappeared 

- all sites i of generation n have renornialized transverse fields h^{i) that can be reduced with respect to their initial 
random fields h™^{i) by a factor r„(i) < 1 (sec Eq. lASp 

h'^(i) ^ /.^(^)r„(^) (2) 

and are connected to the formal external spin via some renormalized coupling J^'^^*(i). These two renormalized 
variables {h^{i), J^'^^'^{i)) take into account the formation of quantum clusters within the decimated generations 
m < n — 1 and the renormalization of couplings with the formal external spin. However the coupling of a site i of 
generation n with its ancestor of generation {n + 1) has not been able to change yet and has still the value J of the 
initial model. 

To renormalize the generation n that is at the boundary, we impose that all sites of generation m > n -\- 1 and all 
external couplings are not decimable. Then the only variables that can be decimated are the renormalized transverse 
fields h^'{i) of the sites of generation n, and their coupling J to their ancestor. Then the renormalized properties 
of a boundary site of generation (n + 1) can be computed from the properties of K independent boundary sites of 
generation n as we now describe. 



B. Properties of a site of generation {n + 1) from the properties of K independent sites of generation n 

We consider a site of generation (n + 1) with an initial transverse field h™+i- It is connected via the initial 
coupling J to K independent boundary sites i = 1,2.., K of generation n, that have renormalized transverse fields 
(*) = ^rr*(0'^n(*) rcnormalizcd external couplings J^'^^*(i) 

Since the site of generation (n + 1) is declared to be undecimable for the moment, we may first renormalize 
independently each branch i — 1,2.., K 



1. RG rule for each branch i = 1,2, .., K 

For each branch i = f , 2, ..K , we apply the Strong Disorder RG rules as follows (see Appendix f^ . We consider the 
maximum between the renormalized transverse field h^{i) — /i™'(i)r„(i) and the coupling J to the site of generation 
(n + 1) 

(a) if hf^{i) ~ h'-^^ {i)rn{i) > J : we decimate the site i (see case (i) in Appendix [A| so that the site of generation 
(n + 1) will get no contribution from branch i to renormalize its transverse field 

(*) =1 (3) 

whereas the contribution of branch i to the external coupling will be (see Eq A2) 

(b) if h!^{i) = h'-^'^ {i)rn{i) < J : we decimate the corresponding bond (see case (ii) in Appendix [A| . i.e. the site i 
of generation n is merged with its ancestor of generation {n + 1). The contribution of site i to the renormalization of 
the transverse field of generation {n + 1) corresponds to the reducing factor (see Eq. A3) 

r-ri^)^^^M^)^^^ (5) 
whereas the contribution of branch i to the external coupling will be 

Jr"'(*) = J„^^^*(^) (6) 
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2. Taking into account the global effect of the RG of the K branches 



We now have to take into account the global effect of the K independent branches to obtain the renormalized 
properties of the site of generation {n + 1) 



K 

Tn+l =n 
i=l 
K 



C. RG equation for the joint probability distribution Pn{J'^^^,r) 



The above rules may be summarized by the following equation for the joint probability distribution P„(J'^^*, r) of 
the renormalized external coupling J^^* = Jn^^*^{i) and the reducing factor r ~ h^{i)/h""{i) of the renormalized 
transverse field for a site i of the boundary generation n 



i—1 



K 

n 

4=1 



0(CV, > J)J(rr'" - 1)5{J^ 



J 



-) + eih^^n < J)5{rr 



K 



K 



S J" 



(5 I r - J]^ r;- 

\ i=l 



(8) 



The first line corresponds to the drawing of K independent sites of generation n with their renormalized variables of 
joint probability Pn{Ji'^* ,ri) and to the drawing of K initial transverse fields /i™' with the law 7r(/i™*). The second 
lign takes into account the renormalization within each branch (j) that give the contribution ( J"^™, r"'^™). Finally 
the third line takes into account the global effect of the K branches. The initial condition is simply 



(9) 



For a finite Cayley Tree containing N generations, the final state after {N — 1) RG steps corresponds to the single 
Center site which will have a renormalized transverse field given by a reducing factor reenter and a renormalized 
coupling J center givcu by taking into account [K +1) independent branches of law P]\f-i{J^'^* ,r) obtained from Eq. [8] 



K+l 



Tjcenter / rext 

\'J center center 

K+l 

n 

i=l 

d(j: 



K+l 



) = n / djr'dnPn{jr\n) I dh^r<^n j dj^'^dri 
i=i ■' ■' 



jnew re 

K+l 



J Lini 

— ^) + 9{hru < mrr-" - n^)5{jr'" - jr) 



ext 



^ I "^cente 



n 



(10) 



D. Pure model without disorder 



To get some idea of the meaning of the RG Eq. [8l let us first consider the pure quantum Ising model with uniform 
transverse field Hq and uniform coupling J. 

The boundary RG procedure described above gives the following results : 

(a) for ho < J, only bonds are decimated, and the whole tree is a single ferromagnetic clusters which is connected 
to the external via the renormalized coupling that grows exponentially in n 

jext ^ J^r^jext (^^^ 
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(b) for ho > J, only sites are decimated so that no quantum clusters are formed, and the transverse fields are not 
renormalized. The external coupling after n generations behaves as 

Jn' =Jn-\(K^)^J^^'(K^) (12) 



So even if each branch yields a reducing factor J /ho < 1, the branching ratio K can overcompensate this one- 
dimensional decay. In this analysis, the pure transition thus takes place at 

jpure ^ ^ ^^3) 

K 

which coincides with the transition found within the 'Cavity Mean-Field' approach (see the discussion on various 
approximation in fx^)- 

The disordered phase ho > KJo then correspond to the exponential decay 

je.t ^ je.t^-^ ^.^^ ^^^^^ ^ 1 OC (Je - J)-^-'- (14) 

In-f 

with the correlation length exponent 

Impure — 1 (15) 

The ordered phase ho < KJo correspond to the exponential growth J^^* = Jq^* ( 7" ) " should stress here that the 

diverging correlation length ^p^r-e defined from J^^* should not be confused with the correlation length ^2 governing 
the decay of the two-point correlation function between two points at distance n 

C(n) cx = e^ft (16) 

The comparison with Eq. [12] shows that the two correlation lengths are related via 

^ ^ InK (17) 



coming from the exponential number if" of points at generation n. In particular at criticality when ^^ure diverges 
£,pure — > +00, the correlation length ^2 takes the finite value ^g*"*** = — !^ agreement with the studies |2TI. [23|. 

Note that for the pure case, the present method based on second order perturbation theory in J/ho will be more 
and more justified in the critical region in the limit of large branching ratio K since the transition takes place at 
smaller and smaller value of the ratio Jc/ho ^ 1/K ^ 0. 



E. Relevance of small disorder around the pure transition 



In finite dimension below the upper critical dimension d < d^, the relevance of small disorder at the transition of 
the pure model is determined by the Harris criterion [23] or the inequality v >2/d oi Chaves et al [2411 . Above the 
upper-critical dimension d > d^ oi the pure model, the Harris criterion is more involved pol [25l . l26l| since there is 
no hyperscaling and there are two distinct correlation length exponents. For the special case of the Cayley tree that 
corresponds formally to d — -t-00, the naive adaptation of the Harris criterion yields that small disorder is always 
irrelevant at the pure transition for the following reason : the correlated volume ^'^ in finite dimension d becomes 
V — on the Cayley tree of branching ratio K\ so the Central-Limit disorder fluctuations which are of order 1/ \/^ 
in finite dimension d becomes exponentially small of order XjsJ on the Cayley tree. So these disorder fluctuations 
are always much smaller than the power-law distance to criticality of the pure model (Jc — J) oc ^-'^l^v^^^^ (Eq. rH| . 

Note however that for quantum disordered models, these arguments based on the weak-disorder Harris criterion 
have been questioned [1, Q since rare regions can play an essential role since the disorder is actually 'infinitely' 
correlated along the time-direction. The effects of rare regions are discussed in detail in the review [23|, with the 
conclusion that the important parameter is the effective dimensionality d^^i of rare regions (in our present quantum 
model where the disorder is actually 'infinitely' correlated along the time-direction, the dimensionality of rare regions 
is rf_R_R = 1 This dimensionality dij^ of rare regions should be compared to the lower critical dimension df"''''^ 
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sufficient to obtain ordering (in our present quantum model where the order is a ferromagnetic magnetization, the 
lower critical dimension is df"'"'' = !)■ The random quantum Ising model thus corresponds to the case d^jj = df"'"'', 
i.e. to the so-called class B in the classification described in [23|, where rare regions dominate the critical behavior 
and can induce an unusual activated scaling. (The conventional power-law scaling is expected to hold for the so-called 
class A corresponding to dun < df , where rare regions cannot undergo the phase transition by themselves, whereas 
the so-called class C corresponds to the case duR > df"^''^ where rare regions can order by themselves at different 
values of the order parameter). 

F. Qualitative description of the possible scenarios for the random critical point 

Before entering the technical solutions of the RG Equation [8] in various phases and regimes, let us now briefly 
summarize the two possible scenarios that we will find for the transition of the disordered model on the Cayley tree 
of branching ratio K > 1 : 

(a) either the random quantum transition takes place in the region where quantum ferromagnetic clusters do not 
exist (as in the pure case described in section IIIDI above) or remain finite. Then the ferromagnetic coupling to the 
boundary will be found to behave in the disordered phase as J^^* ~ e-"/5typ where the correlation length S^typ diverges 
with the typical correlation length exponent Vtyp — 1 which coincides with the pure exponent of Eq. 1151 This case 
can be thus understood as a 'weak-disorder' case where the critical exponents of the pure transition are not changed 
by the disorder, in agreement with the naive adaptation of the Harris criterion discussed in section fll El above. 

(b) or the quantum transition takes place at the point where an extensive quantum ferromagnetic cluster appears, 
and is thus completely different from the pure transition on the tree. We will obtain below that the renormalized 
transverse field scales as ln/i/{(n) cx for n generations, where the distance n* diverges with a correlation 
length exponent of order i'^ ~ 0.75. 

III. BOUNDARY STRONG DISORDER RG FOR THE ONE-DIMENSIONAL CHAIN K = 1 

Before studying the case of real trees having K > 1, we consider in the present section the case of the one-dimensional 
chain corresponding to = 1 to check the validity of the Boundary Strong Disorder RG procedure by comparison 
with the exact solution of the usual full RG equation [l[ . 

A. Choice of Binary disorder 

To simplify the technical details, we now focus on the binary distribution for the initial transverse fields 

7r(/i"") = pid{h"" ~hi)+ P2d{h'"' - h^) with pi+P2 = l (18) 
and we also fix the ferromagnetic coupling to the value J = \/^i^2 i.e. 

hi J hi 



J h2 y h2 

so that the control parameter of the transition is the fraction p2 = I — pi oi Eq. [TH] The exact criterion for the critical 
couphng Jj, [11] 

In Jc = In /i = pi In hi + p2 In /12 (20) 
leads here with Eq. [19] to the critical point 

pr*' = i-pr"-^ (21) 

Starting from the initial condition of Eq [9l it is easy to see that the only renormalized values that can appear are 
the discrete values 

hi 



rfe = ( y ) with k = 0,1,2,.. (22) 



and 



J^^* = Jo^^* ( % ) with m = 0, 1, 2, .. 



J 



i.e. the joint probability distribution P„(J'^^*,r) takes the form of a double sum of delta functions 



+ 00 +00 



jext _ jext 



hi 



hi 



with the normalization 



/»! +00 +00 



m=0 fc=0 

and the initial condition (Eq[5]) 

Cn=Q{m, k) = 5mfi5kfl 

Plugging the form of Eq. [24] into the RG Equation of Eq. [8] yields that the double generation function 

+ 00 +00 



m=0 fe=0 



satisfies the linear recurrence 



c„+i(2/,z) 



piz - 



P2 



Cn{y,z) +P2 



Cn(y,0) 



with the initial condition (Eq 



B. Solution for arbitrary finite-size n 

Let us introduce the generating function with respect to the size n 

+00 

G{y,z,X) = ^A"c„(y,z) 
Summing Eq [28] over n after multiplying by A", yields with the initial condition of Eg [29] 



+ 00 

^ A"c„+i(?/,z) =- 



ji=0 



+ 00 



J2x'^c„,iy,z) 



jn—1 



G{y,z,X) [ 



Piz 



P2 



= -[G(y,z,A)-l] 



P2Giy,0,X) 



i.e. 0(2/, z, A) satisfies the linear equation 



of kernel 



1-A(piz + ^) G{y,z,X) + Xp2G{y,0,\) 

z 



N{z,X) = l- X{piz + ^) 



1 = 



The two roots of this kernel reads 



Z±(A) 



2Xpi 



The sohition which is regular for A is Z (A). Replacing z hy Z (A) in Eq. [32] yields 

1 



Ap2G(y,0,A) 



z-ix) 



y 



-1 = 



G(y,0,A) 



Z-(A) 



Xp2{l-yZ-{X)) 
We may now plug this into Eq|32]to obtain the full solution 



G{y,z,X) 



l-Ap2G(y,0,A) [i^] 



{l-yZ-(\))z 



[l-A(piz+p2|)] i^[^-Z-{X)][z-Z+(X)] {Xpi)[{l-yZ-{X)){Z+{X))-z] 
The expansion in y and z yields the generating function of the joint distribution c„(m, k) 



(m,fc) ^ ^ A"c„(m,fc) = — [Z-(A)]" (^^^ 



n=0 



The joint distribution c„(m, fc) for any finite-size n can be thus obtained by inversion 

dX 1 



dX 



2i7rA"+i 2i7rA"+i Xpi 



1-^/1- 4A^ 



PlP2 



2Api 



2Api 



1 + a/1 - 4A2piP2 



fe+i 



dA 1 
2i7rA"+i A^ 

1 



1 - Vl - 4A2pip2 



2Api 



1 - vi-4AVp2 



2Ap2 



fc+i 



dA 



7c„ 2i7rA«+(™+i)+('=+i) 



1 - - 4A2pip2 



where Cq is a circle around the origin A = in the complex plane. 



C. Asymptotic distribution of renormalized transverse fields as n — ^ +00 

Let us now consider the probability distribution Pn{r) of the variable r alone 

p„(r)= / dr-'p^{r-\r) 



In terms of the representation of Eq. [211 this corresponds to the probability distribution of the integer k alone 

+00 



c„(fc) = ^ c„(m, k) 



m=0 



From Eq. [501 and [371 its generating function reads 

+00 +00 



Y^X"Y.z''cnik)=G{y = l,z,X) = 



n=0 fc=0 



iXp,)[l~Z-{X)] [Z+iX))-z] 
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The series expansion in z yields 

+00 



n=0 



(Api)[l-Z-(A)] [Z+(A)r 



fe+i 



2Api - 1 + Vl - 4A2 



1 - y/l - 4A2p ip2 



-1 fc+i 



so that after inversion, one obtains 



Cn{k) 



dX 



Co 



2Api - 1 + v/l~4A2^ 



1P2 



1 - y/l - 4A2piP2 

2Ap2 



The complex function contains two cuts ] — cx), — Ac] and [Ac, +00 [ on the real axis with 

1 1 



Ac 



In addition, there exists a simple pole at 



> 1 



(43) 



(44) 



(45) 



Ap = 1 ii pi = 1 - p2 < 



(46) 



In this region Pi = 1 — P2 < 5, one thus obtains after deformation of the contour in the complex plane that Cn{k) 
converges towards a finite limit Coo{k) given by the residue of the pole at Ap = 1 yielding 



disordered 



{k) 



n— >-+oo 



disordered 



(k) 



pi 



l-2pi 



if Pi< 2 



(47) 



In the complementary region pi = l—p2 > ^, the variable k does not remain finite as n — > +00, and its asymptotic 
behavior will be governed by the singularities near the branching points (±Ac) of the two cuts of Eq. However, 
to see more directly the appropriate scaling limit, it is more convenient to treat the large integers n and k as real 
variables, and to replace finite sums by integrals, so that the generating function in A < 1 becomes a Laplace transform 
in the variable s = - In A > (Eq H5)) 



+00 



dne '*"c„(fc) 



+00 



E^"c„(fc) 



.n=0 



X=e- 
2 



2pie-^ - 1 + ^1 - 4pip2e-2' 



1 - v^l - 4pip2e 
2p2e-^ 



fc+i 



(48) 



To study the asymptotic behavior of c„(fc) as n — +00, we have to determine the leading behavior as s — >■ of the 
right- handside. 

Let us first consider the critical case pf^*^ — p^''*** — where Eq l48l becomes 



-f CO 



j_ —sn criti 

dne c„ 



{k) 



oc -\/ -e 
s^Q V s 



i.e. after inversion of this Laplace Transform, one obtains that Cn{k) is a half Gaussian defined for k > 



C (/c) cx e(k>0)\j—e-^ 



(49) 



(50) 
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Let us now consider the ordered phase pi = 1 — p2 > 57 where Eg 1481 becomes in the Umit s — ^ 

dne-'''C^'''''"^{k) cx -e L <="-^> <^"-^>' J (51) 

s^o(2pi-l) 

i.e. after Laplace inversion, one obtains that c„(fc) is asymptoticaUy a Gaussian distribution 

1 [fc-(2pi-l),.l2 

cr^-'^(fc) cx ^- e ^Hl^ (52) 



D. Asymptotic distribution of renormalized couplings J*^"^' as n — >■ +cx3 
Let us now consider the probabihty distribution Q„(J'^^*) of the variable J*^^* alone 

In terms of the representation of Eq. I24| this corresponds to the probabihty distribution of the integer m alone 

+00 

qn{m) = y^^qnim,k) 

k=0 

From Eq. [30] and [37l its generating function reads 

+ C30 +C30 

E ^" E = Giy, ^ - 1, A) - 

n— m— 

The series expansion in y yields 

+00 

^A"9„(m) = 



1 



iXp,)[l-yZ-{X)] [Z+iX))^l] 



n=0 



(Api) [Z+{X)) 1] 



1 + v/1 - 4A2pip2 - 2Xpi 



1 - yi - 4A2pip2 
2Api 



(53) 



(54) 



(55) 



(56) 



so that after inversion, one obtains 



qn{m) 



dX 



Co «^A"+1 



1 + ^1- 4A2pip2 - 2Api 



1 - y/l - 4A2piP2 

2Xpi 



The complex function contains two cuts ] — 00,— Ac] and [Ac,+oo[ on the real axis with Eq[ 
exists a simple pole at 



A, 



if pi = 1 - P2 > 



(57) 

In addition, there 
(58) 



In this region pi = 1 — p2 > 5, one thus obtains after deformation of the contour in the complex plane that qn{m) 
converges towards a finite limit Qoo (m) given by the residue of the pole at Ap = 1 yielding 



ordered 



(m) 



CX 
n— >+C30 



ordered 



(to) 



2pi-l fl-pi 



Pi 



Pi 



if Pi > 



(59) 



In the complementary region pi = l—p2 < ^, the variable to does not remain finite as n — >■ +00, and its asymptotic 
behavior will be governed by the singularities near the branching points (±Ac) of the two cuts of Eq. [57l However, 
to see more directly the appropriate scaling limit, it is more convenient to treat the large integers n and to as real 
variables, and to replace finite sums by integrals, so that the generating function in A < 1 becomes a Laplace transform 
in the variable s = — In A > (Eg [56)1 



+CX) 



dne ^"qn{m) 



+00 



.ra=0 



J \=e- 

2 



1 + y/1 - 4piP2e"2^ - 2pie- 



1 - y/l ~ 4pip2e~2" 
2pie-^ 



(60) 
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To study the asymptotic behavior of (7„(m) as n — )■ +00, we have to determine the leading behavior as s — > of the 
right-handside . 

Let us first consider the critical case Pi = P2 — where Eq |60] becomes 



+00 



[1 + Vl~' 



i-VT 



oc W-e-™^ (61) 

s-i-O V S 

i.e. after inversion of this Laplace Transform, one obtains that 9,1 (m) is a half Gaussian defined for k > 

qn"'\m) cx 6l(m > 0)W— e"^ (62) 

n-i.+oo V nn 

Let us now consider the disordered region pi = 1 — p2 < ^ , where Eq |60] becomes in the limit s — > 

+ 00 -, \ s . 2pi(l^pi).^ 1 

dne-'*"g,f^°'-'^'='''='^(m) a L <^-^"' (^-^^i)" J (63) 

s->o (1 - 2pi) 

i.e. after Laplace inversion, one obtains that q-ni'm) is asymptotically a Gaussian distribution 

,. . , 1 [m-(l-2pi),il^ 

qf^'""''^<''''"^(m) cx , (64) 

n-i-+oo y'8Trpip2n 

The duality between c*^°'~^'='-<=''(fc) of Eq. ST] and q°'''^'='''=''(fc) of Eq. ESI between cZ'^^''^'^{k) of Eq. El and 
qdisordevedi^^-^ of Eq. El and the identity between <""(fc) of Eq. EOl and q^"*'ik) of Eq. El are in agreement 
with the duality properties of the one-dimensional model itself [ij . So even if the Boundary Strong Disorder proce- 
dure seems to break explicitly the duality of the model by treating differently the couplings and the transverse fields, 
the rules of Eqs |4] and El are sufficiently symmetric to reproduce dual results as it should. 



E. Dynamical exponent z and Griffiths phases 

Li dimension d = 1, the dynamical exponent z in the disordered phase is known to be determined by the exact 
criterion [29|1 



1 - ( T ) (65) 



h 

For the binary distribution of Egs llSi [TOl that we consider, the dynamical exponent reads 



In^ 

Zexact — 7- — (66) 

Infi^ 



It diverges at criticality. The Griffiths phase is the region near criticality where z > 1. Within the exact strong 
disorder RG approach, the dynamical exponent appears as the coefficient of the exponential decay of the logarithm 
of renormalized transverse-field 

Pilnhn) cx e~*l'"''«l (67) 

In h-R-^ — oo 

Within our boundary strong disorder renormalization, we have found the exponential decay of Eq. 27] for the 
variable k — (see Eq. [52]) : this corresponds to the form of Eq. [57] with the coefficient 



_(i^i) 
In^ 

that coincides with the exact value of Eq. [SB 



^ In 



(68) 
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F. Conclusion for the one-dimensional chain K = 1 

These calculations for if = 1 shows that the Boundary Strong Disorder RG procedure is able to capture correctly 
for the one-dimension chain : 

(i) the exact position of the critical point pf^*^ = i (see Eq \TI} 

(ii) the exact critical exponent tp = ^ governing the scaling of renormalized transverse fields and renormalized 
couplings at criticality fEqs [50l and [62|) 

Inh^ cx Inr oc — fc oc —n'^ with ijj = — 

2 

hiJ''^* oc-moc-n'^ with (69) 

(iii) the exact typical correlation length exponent vtyp = 1 governing the divergence of the correlation length in the 
disordered phase (Eq. l64l) 



disordered { ^\ ^ 



In J;^^* cx - / dmmqf'°'''""""'{m) oc 



S,fyp 



^*yp 1 - 2pi °^ (2 '^^^^ ^typ = 1 (70) 

(iv) The dynamical exponent z coincides with the exact value of the criterion of Eq. 1651 

After this important check for the one-dimensional chain corresponding to if = 1, we now focus on real Cayley 
trees having a branching ratio K > I. It is convenient to consider first the simpler cases J > h^ax and J < /imi„ 
before the more complicated case hmin < J < Knax , where [hmin , hmax] represents the support of the distribution 
TT{hi) of random fields. 

IV. TREE IN THE REGION J > hma^, WHERE ONLY BONDS ARE DECIMATED 

When the ferromagnetic coupling J is bigger than the maximal value hmax of the random fields, it turns out that 
only bonds can be decimated and never sites, so that the RG Eq [8] simplifies into 

i.e. all sites are included into a single quantum ferromagnetic cluster. Its coupling to the exterior reads 

Text T^n Text /'70\ 
Jn = Jo ( ' ^) 

and its renormalized transverse-field satisfies the recurrence 



i.e. ln(r„) is the sum of 



lnr„_i(i) +ln-^ 



(73) 



f(n _ 1 

Mr, = K + K^ + ..if" = K— — - (74) 

if — 1 



i.i.d variables ^In —Jt— ) which are all negative here : the asymptotic distribution is thus Gaussian 



P„(lnr„) (X i-_ e '"^^^Jf"' (75) 



with the averaged value 
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and with the variance 

ijini 7^ l-.ini 

al^M^Var{\n^) ^ -—-Var{\n^) (77) 

In conclusion for J > hmax, the tree is extremely ordered since the whole tree is a single quantum ferromagnetic 
cluster, i.e. this region is very far from the transition. 

V. TREE IN THE REGION J < h^in, WHERE ONLY SITES ARE DECIMATED 

When the ferromagnetic coupling J is smaller than the minimal value hmin of the random fields, it turns out that 
only sites are decimated and never bonds, so that the variable r remains always r ~ I and the RG Eq [8] simplifies into 

Pn{r^\r) = s{r - i)p;;-\r-') (78) 

where the distribution of the external coupling J*^^* alone evolves according to 

p;^i(j-*) = n / djr'PniJr) I dhtMhTis ir^' ' E ) (^9) 

i—l \ i—1 * / 

A. Quantitative mapping onto a Directed Polymer on the Cayley tree 

The evolution of Eq. [79] coincides with the evolution with the length L of the partition function Z£^{f3) of a 
Directed Polymer on the Cayley tree at 'inverse-temperature' /? = 1 

jr = Zf^(/3 = l) (80) 

where the Directed Polymer model 

^r(/3) = Ee^p(-/3 E (81) 

RW \ l<n<L J 

contains effective sites random energies given by 

e(n,i(n)) ==ln/ir(i)-lnJ (82) 

B. Reminder on the Directed Polymer on the Cayley tree 

The Derrida-Spohn solution for the Directed Polymer on the Cayley tree [30| yields that the extensive term of the 
partition function of Eq. |81] is given by 



lnZf^(/3) (X -/3/(/3)L + ... if /3 < /3c 

lnZf^(/3) cx -/3/(/3e)L + ... if P > Pc (83) 



in terms of the following function defined on /3 g]0, +oo[ 



(84) 



and where the location /3c of the freezing critical point is given by the condition 

dpf{/3)y=^.^0 (85) 
In the limit /3 -> 0, the normalization J dhP{h) = 1 leads to the expansion 

/(/3) cx -InJ-i^^-oo (86) 
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In the limit /3 — >■ +c», we have to distinguish two cases 

(i) if the distribution P{h) is continuous, one may perform a saddle-point calculation at /imm with the following 
change of variables h — hmini^ + 2;//?) to obtain 



dhh-^P{h) I h„,,^^hJ,^e-^P{h^,r.) "-"^^"""^ (87) 



that leads to the following expansion 



(X — In J + In h 



In 13 \nK + \n h^in + In P(/i,„i„) 



/3-++00 /3 ^ 

The positivity of the leading correction > shows that the positive limit 

/(/3 = +c») = -lnJ + ln/i„„>0 (89) 
is reached from above. Between these two limits, the maximum reached at /3c is thus positive 

/(/3e) > fW = +«))> (90) 
(ii) if the distribution P{h) is the binary distribution of Eq. [18] 



J dhh-^Pih) = pih^^ ^1 



Pi 



(91) 



then the asymptotic behavior reads 



f{f3) oc In/ii-lnJ (92 

/3^+oo fj P Pi 

So the positive limit 

/(/3 = +00) = -In J + ln/ii > (93) 

will be reached from above if Kpi < 1 and from below if Kpi > 1. For the case Kpi < 1, one has thus a finite location 
/3c satisfying Eq[5ni For the case Kpi > 1, there is no freezing, i.e. in the limit of zero temperature /3 — )■ +00, the 
DP is able to find a path containing only hi (since pi is above the percolation transition pi > pP'^'''^° — 1/ K). 



Application to the variable J*^^' of the quantum model 



The solution of Eq. [83] can be directly translated for the external coupling J£^* via the mapping of Eq. [8OI 

Injr* = lnZf^(/?= 1) oc -![1)L+... if 1< /3, 



Injr' = lnZi:'^(/?= 1) cx -J{/3c)L + ... if 1 > 



(94) 



For the case 1 > /3c, the inequality f{/3c) > shows that In J£^* decays with L so that the quantum model can then 
only be in its disordered phase. 

For the case 1 < /3c, we see that the quantum model can be either disordered if /(I) > 0, ordered if /(I) < and 
at criticality if /(I) = 0. This condition for criticality reads using Eq [5H 



1 = KJc J dhh-^P[h) 



The condition 1 < /3c to be in the delocalized phase of the Directed Polymer is equivalent to the condition 

j dhh'^{\nh)P{h) 



< /'(/3 = 1) = In 



K / dhh-'P{h) 



Jdhh-^P{h) 



(95) 



(96) 
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Using the transition condition = f{(3 = 1) = — In [KJc J dhh ^P{h)\ , we may rewrite the condition as 

o</ (/^-i)--hij,+ j^f^f^^.p^f^^ jdhh~ip(^h) ^^^^ 

which is always satisfied if Jc < hmim which is necessary to have the quantitative mapping onto the DP (see Eq. I79p . 



D. 



Conclusion 



The conclusion of this analysis of the region < J < hmin, where there exists a quantitative mapping onto a 
Directed Polymer on the Cayley tree is the following : 

(i) either the transition of the quantum model takes place in this region < Jc < ^mm- Then this can only 
correspond to the delocalized phase of the Directed Polymer, where the external coupling has for asymptotic behavior 
around Jr 



jext 



Z^P{P = 1) oc e-^(i)^ 



KJ j dhh-^P{h) 



In particular, in the disordered phase J < the typical external coupling decays exponentially 



In Jf 



L 

S,typ 



where the correlation length ^jj,p diverges as 



(98) 



(99) 



^typ 



In 



'Jo' 
< J 



OC {Jc - J) with Vfy 



(100) 



Note that this scenario where the ordered phase appeared at Jc < hmin where one decimates only sites and no bonds, 
i.e. where no clusters are formed is completely different from what happens in d = 1 (see section . and is due to 
the branching with ratio if > 1 of the tree. 

(ii) or the quantum model remains disordered in the whole region < J < i.e. J£^* decays exponentially in 



L even for J 



Then the transition will take place in the region h„ 



< J <K 



where one has to decimate 



both sites and links. This is the case studied in section IVT 



E. Special case of the binary distribution 

For the binary distribution of Eq. [18] discussed around Eq. [93l various cases can appear : 

(a) For the case pi > there is no freezing at finite /3c for the corresponding Directed Polymer model, so that the 
behavior of the external coupling is given by Eq. IHHlin the whole region < J < /imin = hi 



jext 
J J. 



Z^P{P = 1) OC e-/(i)^ 



(101) 



where 



K J dhh^^P{h) 



El. 

hi 



P2_ 

h2 



< 



hi 



(Kpi) 



< hi 



(102) 



i.e. the quantum transition always takes place at Jc < hmin — hi, where there exists the quantitative mapping onto 
a Directed Polymer model in its delocalized phase. 

(b) For the case pi < there exists a freezing transition for the associated Directed Polymer model at some finite 
/3c, and one has to discuss whether /3c > 1 or /3c < 1. 

According to the discussion of section III Dl a pure system with the uniform random field ho taking the value hmax 
would have a pure transition at J^^^"^ {hmax) = so one expects a lower critical coupling for the disordered case 

when some random field are smaller Jc < J c^^'^ {hmax) = '^"i^" ■ So if one chooses ^"t?" < hmin, one is sure to have 
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J < hmin, i-e. the transition will takes place in the region where only sites are decimated (section |V| and where there 
exists the quantitative mapping onto the Directed Polymer described above. 

On the other hand, if one chooses a binary distribution {hi < /12) with pi — I ~ p2 ^ 0, then the transition will be 
near Jc ~ J^"'''^(/i2) = so if one chooses /ii < the transition should take place in the region hi < J < h2 where 
one needs to take into account the formation of quantum ferromagnetic clusters, that we discuss in the following 
section. 

VI. TREE IN THE REGION < J < h^a^, WHERE BOTH SITES AND LINKS ARE DECIMATED 

In this section, we consider the remaining case hmin < J < hmax where both sites and links are decimated. 

A. Evolution of the distribution of the variable r 

Let us first focus on the distribution of the renormalized transverse fields alone, i.e. on the distribution Pn{i') 

P„{r) = J dr-'P„ {r-\r) (103) 
that evolves according to the recurrence (after integrating Eq[5]over J^^*) 

P„+i(r) - n / dr,P„(rO / dr^^ [ dJC'irih^^) 



9{h^rn > J)^(rr" - 1) + OihY^'n < J)S{rr'" - n^) 



with the initial condition 



(104) 



Pn=o{r) = 5{r - 1) (105) 



B. Special case of the binary distribution 

To simplify the technical details, we now focus again on the binary distribution of Eq. [18] and we fix /ii < J = 
\/hih2 < /i2 so that Eq [19] holds. The control parameter of the transition is the fraction pi — 1 ~ p2 oi Eq. [18] 

Starting from the initial condition of Eq 11051 the only values that can appear are of the form of Eq. [52] i-e. the 
probability distribution P„ (r) takes the form of a sum of delta functions 



+00 



P„(r) =^c„(fc)J 



k=0 



, \ k' 

hi 



(106) 



with the normalization 



1= / drPn{r)^Y.''^('') (10^) 
•^0 k=0 



and the initial condition (Eg 1105 



c„=o(fc) = 4,0 (108) 
Plugging the form of Eq. 11061 into the RG Equation of Eq. 11041 yields that the generation function 

+00 

c„(z) = ^z'=c„(fc) (109) 

A:=0 
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satisfies the recurrence 

K 



C„+i(z) = ( C„(z) piZ+— +P2C„(0) 



1-i 

z 



(110) 



with the initial condition (eq[ 

Cn=0{z) = 1 (111) 

In one dimension _fC = 1, where the corresponding equation was Unear and could be solved for any n via the 
introduction of the time generation function (Eg I31|) . Here for A' > 1, the equation is not linear anymore, and we 
do not know how to solve exactly Eq. IllOl for any finite size n. So the goal will be here to determine the possible 
asymptotic behaviors of the distribution c„(fc) as n — > +oo, as a function of the parameter pi = 1 — P2 g]0, 1[. 

C. First discussion on the possible behaviors as n ^ +cx3 

It is clear that the asymptotic behaviors as n ^ +oo will depend on the behavior of Cn{k) near the origin k = 0,1,2,.. 
and on the behaviors of the first moments. It is thus convenient to write the corresponding recurrences. 

1. Recurrences for Cn{k) near the origin k — 0,1,2 

From Eq. 11101 one obtains that the recurrences for c„(fc) for the first values of /c = 0, 1, 2 reads 

c„+i(0) = [P2(c„(0) + c„(l))]^ (112) 

C„+i(l) - K b2(c„(0) + Cn{l))f'' [P1C„(0) +P2C„(2)] 

C„+i(2) = K b2(c„(0) + C„(l))]^-' blC„(l) +P2C„(3)] + '^^^^^ [P2(C„(0) + C„(l))]^-' [piC„(0) +P2C„(2)]' 

In the limit n — > +oo, there are only two possibilities : 
(i) either Cn(fc) converge towards finite values Coo(fc) 

c„(fc^0,l,2,..) cx Coo(fc = 0,l,2,..) >0 (113) 



that satisfy Eqs 11121 as stationary equations. 

(ii) or Cn{k) for finite values fc = 0, 1,2, .. converge towards zero as n +oo : then from Eq. I112| the only way 
they can decay to zero is exponentially with respect to i^T" 



c„(fc = 0,l,2,..) oc e-(^^*)^" (114) 



2. Recurrences for the first moment < k >„ 
Using the series expansion in the moments 

Cn{z = e-") = J2 e^^'cfc = E ( 1 - '^^^ + ^ + ^(9') J Cfe = 1 - g < A; >„ +|- < fc2 ^^^^3) (^^5) 

k=0 k=Q ^ ^ 

one may derive from Eq. IllOl the recurrences for the first moment 

< k >„+i^ if [< k >„ +(2pi - 1) +p2C„(fc = 0)] (116) 
In the limit n — ^ +00, there are only two possibilities : 

(a) either < k >„ converge towards a finite value < k >oo that satisfies Eq 11161 as stationary equation, where 
Coo(fc = 0) has to remain finite 

< k >oo= K[<k>oo +(2pi - 1) + P2Coo(fc = 0)] (117) 
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For K = 1, this fixes the value Coo{k = 0) = (1 — 2pi)/p2 found previously (Eg |47|) whereas for K > 1, it yields 

<k>^^ -^[{l-2pi)-p2Coo{k^0)] (118) 

which should be positive, since A: is a positive random variable, so we have the bound 

< Cooik = 0) < ~ ^^'^ implying p, < ^ (119) 
P2 2 

In fact, since k is an integer random variable, we have in fact the stronger constraint 

+00 +00 

< k >oo= ^ fcCoo(fc) > ^Coo(fc) = 1 - Coo(fc = 0) (120) 
k=l k=l 

which yields with Eq. 11181 

(l-2i^pi)-(l-i^pi)coo(fc = 0)>0 (121) 

For instance ioi K = 2, this yields 

cl^=^\k - 0) < implying < ] (122) 

1 - 2pi 4 

(b) or < fc >n diverges as n +oo : then Eq. II 161 can be approximated by < fc >„+i~ K < k >„ so the only way 
it can diverge is like 

<k>n oc if" (123) 

n— ^+oo 

3. Recurrences for the variance o"^ =< >„ — < fc >„ 

From the recurrence concerning the second moment 

< fc2 >„^,= K[< fc2 >„ +(i;c - 1) < fc >2 +2 [K{2pi -1) + {K- l)p2c„(fc = 0)] < fc >„ 

+if (2pi - 1)2 + 4pip2 + P2 (2X(2pi - 1) + 1 - 4pi) c„(fc = 0) + (i^ - l)p^c2 (fc = 0)] (124) 

one obtains that the variance 

al =< fc2 >„ - < fc >2 (125) 

satisfies the recurrence 

<jI+, = K [gI - 2p2C„(fc = 0) < fc >„ +4pip2 + (1 - 4pi)p2C„(fc = 0) - (fc ^ o)] (126) 
In the limit n — )■ +cx), there are only two possibilities : 

(a) either cr^ converge towards a finite value cr^ should be stable by the recurrence of Eq. 11261 

'^lo =i^['T^-2p2Coo(A: = 0)<fc>oo+W2 + (l-4pi)p2Coo(fc = 0)-p2c^(fc^0)] (127) 
i.e. using Eq. I118| this yields 

K 



2 



{K- 



--^P2 [-P2{K + l)cL(fc = 0) + (i^ + 1 - 4pi)coo(fc = 0) - 4pi{K ~ 1)] (128) 



The positivity of the variance implies that the discriminant of the right-handside second order polynomial should be 
positive, and that Cao(k = 0) S]c_, c+[ where c± are the two real roots. 
For instance for K = 2^ Eq. 11281 becomes 



-2p2h3p2C^(fc = 0) + (3-4pi)coo(fc = 0)-4piJ (129) 
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with the discriminant 

S = 9- 72pi + &Ap\ (130) 
Since < pi < 1/2, this discriminant is positive only if 

Pi < pT" - 3(3 - ^/5) ^ 0^^43237 ^^3^^ 
16 

(b) or CT^ diverges as n — > +oo : then from Eq. 11261 the only way it can diverge is like 

al cx (132) 

4- Summary 

From the possible behaviors of c„(A; = 0) and of the two first moments as n — )■ oo, we have seen that either a finite 
c„(fc) exists with no rescaling in k, but this can only be the case for sufficiently small pi, or the variable k flows 
towards infinity with the first moment and the variance growing as if", as we already found for J > h^ax in section 
IIVI Let us first discuss the case pi > ^ where the solution found in section HVl can be simply extended. 

D. Asymptotic rescaled distribution Coo(fc) for pi > i 



For pi = 1, one has J > h^ax — hi and one should recover the solution of section HVl It is clear that this type of 
solution can survive in the whole region pi > ^ for the following reasons. Let us assume that c„(0) is negligible in Eq. 
IllOl : then, taking into account the initial condition of Eq. Illli one has the simple solution in terms of M„ of Eq. [73] 



c„(z) = [c„_i(z) piz + — ] = (co(z))^ piz+— =[piz + — ] ^Y^Clj^[pizf — (133) 
i.e. the coefficient of order k simply reads fEg 11091) 



J=0 



Mn+k Mn+k M^-k 

Cn{k)^C^,; Pi = (134) 
For n — > +00, we may use the Stirling formula to obtain the asymptotic Gaussian formula 

1 [fc-(2pi-l)M„l^ 

that coincides with Eq. [55] for the one dimensional chain K = 1 where Af„ = n. This solution is consistent as long as 

c„(fc = 0) (X e-^^-^Hi^ (136) 

remains exponentially negligible in Mn with respect to the center of the distribution localed around k ~ (2pi — l)Af„ 
i.e. for pi > i. 

E. Existence of a finite asymptotic distribution Cooik) as n — >^ 00 for small enough pi 

If a finite asymptotic distribution Coo(fc) exists as n — ?> +00 without any rescaling in /c, it should satisfy the stationary 
version of the recurrence of Eq. IllOl 

z'^C^iz) ^ {c^{z)[piz'^ +P2]+P2Coo{Q)[z-l])^ (137) 

Note that for pi = 1 — p2 = 0, we know that we have the trivial solution 

ct^-^){z)^l i.e. c(2i=o)(fc)=4,o (138) 
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A stationary solution Coo(fc) is thus likely to exist for small enough pi = I — P2- 

Let us now focus on the case of branching ratio K — 2, where this equation is simply quadratic 

[piZ^ + P2]^ Cl{z) + (2p2Coo(0) [piZ^+P2] [Z~l]~ Z^) Coo (z) + P^C^ (0) [z - if = 

and involves the discriminant 

A(z) ^ z^^pIc^{Q)R{z) 



R{z) = 



l + P^z' 

P2 



[1-z] 



4picoo(0) 



The two possibles solutions read 



(2p2Coo(0) [piz^ +P2] [l-z] + z2) ± yA(I) 



2[piz^+P2r 

At 2; = 0, both solutions coincide with Coo(O). At z = 1, the value A(z = 1) = 1 yields 

= 1) = 

So the solution with the proper normalization 1 = £00(2: = 1) is the solution with the (+) sign 

Riz) 



Coo(z)=C+(z) =Coo(0)- 



1 + I^z^ 

P2 



Coo(O)- 



(139) 



(140) 



(141) 



(142) 



Riz) 



-R{z) 



(143) 



Since R{z) of Eq. 11401 is polynomial of third degree having R{z = 0) = 1, these expressions shows that Eg 11431 can be 
expanded in power-series of z, and the coefficients Coo(fc) can be computed via the following contour integral in the 
complex plane 



Coo (^) 



dz 



:yCoo(2:) 



(144) 



Now we have to impose that these coefficients c„(/c) that represent probabilities are all real and positive. So we have 
to discuss the cuts that appear in ^jR{z) i.e. the three roots Zi of the cubic polynomial R{z) of Eq. 11401 that can be 
written as 



R(z) = az^ + bz'^ + cz + d ^ — TT(zi - z) 

P2 



(145) 



with parameters 



d = 1 

c = -1 

b =^P' 



1 



P2 4ppoo(0) 

Pi 

a — 

P2 

To determine the three roots Zi of Eq. 11451 it is convenient to perform the following change of variables 



(146) 



Z = X = X 

3a 



1 



4piP2Coo(0) 



so that the new variable x satisfies the standard form of cubic equation 

x^ +Px + Q = 



(147) 



(148) 



21 



of parameters 

52 _ [l+8piP2Coo(0) + 16pip2(4pi-3)c2 



P ' 



d b 

Q = - H 

a 27a 



26^ 9c 
a 



[1 + 12pip2Coo(0) + 24pip2(5pi - 3)c^(0) + GiplpU^ - 8pi)c] 



,3 



864p3p3c3^(0) 



(149) 



The discriminant of this cubic equation reads 



^3 ^ 



A, -0^+27^' 432ptp2£3^(0) 
N3 = 1 + P2(13pi - l)coo(O) + I6plpi{8pi - 5)c^(0) + 64pip3c^(0) (150) 

The numerator A'3 is itself a third-degree polynomial in the variable c^(0) ; again it is convenient to perform the 
change of variables 

£oo(0)=i+|--^ (151) 
3 4p2 

to rewrite the numerator in the reduced form 

N3 ^ Gipiplit^ + pt + q) (152) 



of parameters 



P 



3 + 64pi(l-4pi) 
192pip2 

63 - 280pi + 3584pf - 4096pj 
6912]5ip2 



(153) 



The corresponding discriminant reads 



^3 -q^ + ±p^--^'-''''h (154) 
^ ^ 27^ 1769472pfp4 v y 



The possible solutions Zi of the cubic polynomial R{z) of Eq. 11451 depend on the sign of the polynomial A3 of Eq. 

m\ 

(i) If A3 > 0, there exists a single real root, and two complex roots that are complex conjugates. 

(ii) If A3 < 0, the three roots Zi are real. 

(iii) If A3 = 0, there exists a simple real solution Zg and a double real solution Zd 

Moreover, since the generating function c(z) contains only non-negative coefficients c(fc) > 0, we known from 
Pringsheim theorem that its closest singularity to the origin has to be on the real axis. 

In the limit pi — > where we know the solution (Eq. 11381) . R{z) reduces to the second order polynomial 

i^(Pi=o)(^)^ (2_f)^ (155) 

corresponding to the limiting case where the double real solution is finite, and the simple real root diverges 

(pi=0) _ „ 

z^J"=^^ +00 (156) 

By continuity in pi, we expect that Zs will become finite but large. The double real solution could in principle 
separate, but they cannot become complex (case A3 > 0) as a consequence of Pringsheim theorem, since they are 
closest to the origin than Zg. It turns out that they cannot separate along the real axis either (case A3 < 0), because 
the integrals of Eq. 11441 would have a negative sign. The only remaining possibility is thus that for small enough pi , 
there remains a double real solution Zd corresponding to the condition of the vanishing determinant 



A3 = 



(157) 
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Then the simple real solution Zs and the double real solution Zd read 

sg h 



Zd 



P 3a 
3Q h_ 

^2P ^ 3a 



and R{z) can be rewritten as 



R{z) = ^{z,-z){z-Zdf 

P2 



(158) 



(159) 



i.e. there exists a single cut [z^, +oo[ on the real axis. After the deformation of the contour in the complex plane, the 
coefficients of Eq. 11441 can be written as the real integral 



Co 



dz _ y/p2PlCoo{0) 

2tTTz'^+^''°^^^' ^ TT 



- Zs{x - Zd) 

ax 



x''{p2 +PlX^) 



(160) 



which is positive for Zs > Zd, which should be possible for small enough pi since one starts from the values of Eq. 
11551 In particular, the asymptotic behavior for large k reads 



^ \/p2PlCoo{0)T {^) jZs - Zd) 
fe4^oo 7r(p2+Pl^?)22fe-3/2^3/2 



(161) 



Let us now study more precisely the condition A3 = of Eq. 11571 : this means that Coo (0) has to be one root of 
the cubic polynomial N3 of Eq. 11501 Let us start from the simple case pi = I — P2 — 



^i^^=°^ =1-Coo(0) 



(162) 



that vanishes for Coo(O) — 1, which is indeed the correct solution (Eq 11381) . In the region < pi < the discriminant 
63 of Eq. [154] remains negative, so that the three real roots of Eq. 11521 can be written using the integers m = 0, 1, 2 as 



2\l —— cos 



9 + 2nm 



where the angle < 9 < tt satisfies 
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(63 - 280pi + 3584p? - 4096p?) 



2 V ~P' 

In the limit pi = 1 — p2 ^ 0, we have at first order in pi 



Pi 



P2(3 + 64pi(l-4pi))3 



(163) 



(164) 



P 



64pi 92 
7 77 



0{Pi) 



768pi 3456 
cos 6* = -7^3^1+0(^1^^) 
=l + 7V3fi + 0{pT) 

so that the only root that remains finite as pi — >■ corresponds to the value to = 2 in Eq. 11631 

7 7 



to 



12 



whereas the two other solutions flow towards (±00) 

to 
ti 



:Pi+0{pi) 



(165) 



(166) 



8v^ 24 
-8^ -I 



(167) 
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The 'physical' solution for Coo(O) fEq I15ip has thus for expansion 



Coo(0)-C2=t2 + %-^ = l-2pi + +O(p?) (168) 
The two other unphysical roots of A3 read 

(8p2-3) 1 ^ 1 1/2. 

-^^^Hr^--V^ + \ + 0[A") (169) 
12p2 8^pr 8 

By continuity in pi , we thus expect that Coo (0) = C2 in the whole region < pi < ^ where the three real roots of 
A3 remain distinct. At the border Pi ~ the discriminant ^3 vanishes : there exists a double root Coo(O) — 49/54 
(coalescence of the two solutions C2 et Cq) and a simple root Ci = —16/27 that is negative. For Pi > the 
discriminant 63 becomes positive, i.e. the double root Cq = C2 = 49/54 become a pair of complex conjugated roots 
(this cannot correspond to any physical Coo(O) which has to be real), and the simple root Ci remains real (but this 
cannot correspond to any physical Coo(O) which has to be positive). 

So we have found that a stationary solution Coo(fc) exists in the region < pi < Let us consider more precisely 
the solution at the frontier Pi — where we have found Cqo (0) = 49 /54 : it turns out that the parameters of Eq. 11491 
both vanish at this point P = = Q, i.e. there exists a single triple root to R{z) (merging of the simple solution Zs 
and of the double solution Zd existing in the region < pi < ^) 

zs = zd = 3 (170) 

Then Eq 11601 simplifies into 

4 72 r+°o (x _ 3)3/2 

coo(fc) =-7=/ dx I J (171) 

It is easy to check that it satisfies the value at fc = 

c (k-0) ^ ^ - ^ (172) 

''^"^^ - nV2 J3 (27 + x^r ~ 54 ^^^^^ 



and the normalization condition 



k=0 



In particular, the asymptotic behavior for large k of Eq. I171l reads 



coo(fc) oc V^^^^^.f (174) 

+ piz2)2^^*^-5/2^5/2 



in contrast to Eq. 11611 

It is thus interesting to mention now the singularities that appear for 



We find 



Pi^ -^-<^ with e > (175) 
28 



3 1078 

1 539 

=4 + ^^ + - 

4519603984 . 

d-i = e'^ + ... 

14348907 

36156831872 . 
14348907 ' 

2187^3 
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and 



Coo(0) = C2 

Ci 

Co = 



49 3430 



54 729 
16 6272 



134456 
6561 



27 729 
49 3430 134456 



3/2 



3/2 



(177) 



54 729 6561 

Finally, the parameters P and Q (Eq ll49p and the roots (zs, Zd) fEq ll59p of the cubic equation in z (Eq. I145|) . read 

784 
-—. + ... 

43904 



P 

Q 



.3/2 



= 3 



27 

56 1 



/2 



(178) 



In particular, the difference {zs — Zd) between the simple and the double root presents a square-root singularity in 
y/e. Let us now discuss what happens on the other side of pi = 1/28. 

F. Region Pi ~ -ks + (■ with e > 




100 




40 I 

0.038 



FIG. 2: Region Pi = ^ + e with e > (a) Cn(k = 0) as a function of n for various values of p\ : there is a rapid convergence 
towards a finite value, and then a decay towards zero at some n*(pi) (b) Divergence of the scale n*(pi) as pi — >■ 1/28 = 0.0357: 
the inset displays the same data in the log variables lnn*(pi) as a function of In (pi — ^). We measure a slope of order 
u+ ~ 0.75 (see Eq. [T79l) . 



For Pi — + ^ with e > 0, as shown on Fig. [5] we find numerically from Eq. IllOl that Cn{k) rapidly converges 
towards a quasi-stationary metastable distribution c* (fc) that disappears after a number of generation (e) that 
diverges as e — > as 

n*(pi) oc T—ir- with I/+ ~ 0.75 (179) 

P,^(^)+ (Pi-^) 

This diverging length represents the radial distance near the boundary where the variable k remains finite, i.e. where 
the quantum ferromagnetic clusters remain finite. But then after n*, the variable k flows towards infinity, i.e. an 
extensive ferromagnetic cluster is formed. For a tree of n generations, the renormalized transverse field of this central 
cluster will scale as 



hihii{n) cx —k cx —K' 



(180) 
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This should be compared to the scaling \nh]i{L) oc — in finite dimension d. The corresponding renormalized 

external ferromagnetic coupling will scale as 

J^'^^* ~ if"-"*J^;^* (181) 

where J^^^^ corresponding to the first rt* generations has the behavior of a modified Directed Polymer partition 
function as we explain in the following section. This should be compared to the scaling J£^* cx aL'^^^ in finite 
dimension d > 1 where there exists an underlying classical transition (see {Bl ij and references therein for the properties 
of the classical Ising model on the Cay ley tree). 



G. Modified Directed Polymer model for J'^^* when there exists a finite Coo(fc) 

In the region J < hmin where only sites are decimated and where the asymptotic distribution Coo(fc) reduces to 
the delta function Coo(fc) = 5fe,0i we have seen that the renormalization of the external coupling J^xt corresponds 
quantitatively to a Directed Polymer model described in section |Vl In the region hmin < J < hmax where the 
asymptotic distribution Coo(^) is not a delta function anymore but remains a finite distribution with no rescaling in 
k, the renormalization of the external coupling J^xt still corresponds to a Directed Polymer model, but with a slightly 
different disorder with respect to the case described in section |V] : the effective random energies of Eq. |52] that were 
both positive for the binary case in the region J < hmin, with ci = In ^ with probability pi and 62 — In with 
probability p2, are replaced in the region hmin < J < hmax when there exists a finite Coo(fc) 



with probability p'l = 1 — p2Cco{k ~ 0) 

In — ^ = In \ J— with probability P2 = P2Coo{k — 0) 
J \ hi 



(182) 



The factor Coo(fc — 0) enters because the non-trivial contribution to the renormalization of Jext occurs only if the 
corresponding site gets decimated, i.e. one needs both conditions hi = /12 and k = (see the renormalization rules of 
section IiTbT|) . 

The function of Eq. [M] has thus to be modified into 



fmodiP) 




{1-P2coo{k = 0)) +P2Coo(fc = 0) ( 5^ 



(183) 



Again for 1 > /?c, the quantum model can then only be in its disordered phase, whereas for 1 < /3c, the criticality 
condition fmod{P — 1) for the quantum model of Eq. 1951 becomes 



K 



(1 -P2Cryo{k = 0)) +P2Coo(fc = 0) 



hi 



(184) 



and can only occur within the delocalized phase of the Directed Polymer. 



H. Summary 

The analysis of the region hmin < J < hmax can be summarized as follows : 

(i) either the quantum transition occurs in the region where there exists a finite asymptotic distribution Coo(fc) 
describing finite quantum ferromagnetic clusters : then the transition is analog to the transition that can occur in 
the region J < hmin (see section |V| and is determined by the modified Directed Polymer model described just above. 
The critical behaviors are then governed again by EgslMland Eq llOOl with the critical exponent vtyp = 1- 

(ii) or the quantum model remains disordered as long as there exists a finite asymptotic distribution Coo(fc), i.e. 
up to the critical point pi = 1/28 described above. It becomes ordered only when there appears an infinite quantum 
ferromagnetic cluster for pi — l/28-|-e. Then the critical behaviors are completely different from (i) and are determined 
by the scaling of section IVIFI : the finite-size renormalized transverse field and renormalized external coupling scale 
as Eq. 11801 and II8II where the diverging length n, of Eq. 1 1791 involves the critical exponent z/+ ~ 0.75. 
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VII. CONCLUSION 



In this paper, we have considered the Random Transverse Field Ising model on a Cayley tree. To avoid the difficulties 
of the usual Strong Disorder RG that destroys the tree structure, we have introduced a modified procedure called 
'Boundary Strong Disorder Renormalization' that preserves the tree structure, so that one can write simple recursions 
with respect to the number of generations. We have first checked that this modified procedure allows to recover 
exactly most of the critical exponents for the one-dimensional chain. Then we have studied the RG equations for the 
Cayley tree with a uniform ferromagnetic coupling J and random transverse fields within the support [hmimhmax]- 
We have found the following picture : 

(i) for J > hmax, only bonds are decimated, so that the whole tree is a quantum ferromagnetic cluster 

(ii) for J < hmin, only sites are decimated, so that no quantum ferromagnetic cluster is formed, and the ferromag- 
netic coupling to the boundary corresponds quantitatively to the partition function of a Directed Polymer model. This 
relation with the Directed Polymer has been already obtained for the Cayley tree via some approximations within 
the Quantum Cavity Approach [l7l-[l9j. and for arbitrary networks via simple perturbation deep in the disordered 
phase |20j . However here we also consider the possibility of bond-decimations to build quantum ferromagnetic clusters 
which may become important near the transition. 

(iii) for hmin < J < hmax, both sites and bonds can be decimated. When the quantum ferromagnetic clusters 
remain finite, the physics is similar to (ii), with a quantitative mapping to a modified Directed Polymer model; 
otherwise an extensive quantum ferromagnetic cluster appears. 

We have found that the quantum transition can be of two types : 

(a) either the quantum transition takes place in the region where quantum ferromagnetic clusters do not exist or 
remain finite. Then the ferromagnetic coupling to the boundary behaves in the disordered phase as J^^* ~ e""''^*'"' 
where the correlation length ^typ diverges with the typical correlation length exponent vtyp = 1- 

(b) or the quantum transition takes place at the point where an extensive quantum ferromagnetic cluster appears 
at the center of the tree, at a radial distance from the boundary that diverges with the correlation exponent 
v+ ~ 0.75. In the ordered phase, the finite-size renormalized transverse field and renormalized external coupling scale 
as InhRin) cx -/C"-"* and J^'^^* ~ JC"""* J^^'^^*. 

In a companion paper [s^ , we describe how the idea of Boundary Strong Disorder RG can be adapted in dimension 
d = 2 and we present the corresponding numerical results. 



Appendix A: Reminder on Strong Disorder RG rules on arbitrary lattices 

In this section, we recall the standard Strong Disorder Renormalization for the Random Transverse Field Ising 
Model of Eq. [H 

For the model of Eq. [U the Strong Disorder RG rules are formulated on arbitrary lattices as follows [1, 01 '■ 
(0) Find the maximal value among the transverse fields hi and the ferromagnetic couplings Jjk 

— max [hi, Jjk] (Al) 

i) If Q ~ hi, then the site i is decimated and disappears, while all couples (j, k) of neighbors of i are now linked via 
the renormalized ferromagnetic coupling 

-Jjk - Jjk H (A2j 

ii) If ri = Jij, then the site j is merged with the site i. The new renormalized site i has a reduced renormalized 
transverse field 

h^'^^h.n with r,^^ (A3) 

Jij 

and a bigger magnetic moment 

Mr" = M.+M, (A4) 
This renormalized cluster is connected to other sites via the renormalized couplings 

jrr = J.k + j,k (A5) 
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(iii) return to (0). 
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